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Abstract 

Thermospheric  mass  density  variation  is  the  major  source  of  aerodynamic  drag 
estimation  errors  at  altitudes  below  about  700  km.  This  project  demonstrates  how  the 
current  limit  of  thennospheric  mass  density  predictability  can  be  extended  by 
systematically  integrating  observations  into  a  coupled  thennosphere-ionosphere  first- 
principles  model.  An  ensemble  data  assimilation  procedure,  constructed  with  the  NCAR 
Data  Assimilation  Research  Testbed  (DART)  and  the  NCAR  Thennosphere-ionosphere 
Electrodynamics  General  Circulation  Model  (TIEGCM),  can  take  advantage  of  the  tight 
coupling  between  the  ionosphere  and  thermosphere,  enabling  the  inference  of 
thennospheric  temperature  and  compositions  from  abundant  GPS-based  ionospheric 
observations.  Observing  system  simulation  experiments  suggest  that  thennospheric 
states,  particularly  temperature,  can  be  well  inferred  by  assimilating  electron  density 
profiles  obtained  from  the  COSMIC/FORMOSAT-3  mission  into  the  TIEGCM.  This  in 
turn  leads  to  a  significant  improvement  of  the  neutral  mass  density  forecasting  longer 
than  3  days.  Furthennore,  validation  of  assimilation  analyses  with  independent  CHAMP 
mass  density  observations  confirms  that  the  approach  indeed  improves  the  thennospheric 
mass  density  specification.  Predictability  of  the  ionosphere  can  be  extended  considerably 
by  estimating  neutral  compositions  from  GPS-based  ionospheric  observations  by  the 
coupled  thennosphere-ionosphere  data  assimilation  approach  developed  in  this  project. 


Background 

Even  though  the  Earth’s  upper  atmosphere  density  is  tenuous,  it  is  substantial  enough  to 
exert  significant  drag  on  orbiting  spacecraft  and  debris.  Understanding  this  drag  is 
complicated  by  the  fact  that  this  region  is  subject  to  highly  variable  external  drivers  from 
above  and  below.  Also  its  ionized  constituents  affect  telecommunication  and  navigation, 
motivating  numerous  observational  and  modeling  efforts  since  the  dawn  of  space 
exploration.  Because  of  a  lack  of  global  observations  of  thennospheric  parameters,  our 
understanding  of  the  density  variability  rests  largely  on  the  averaging  of  numerous 
observational  samples,  limiting  our  modeling  capability  to  the  realm  of  climatology.  To 
overcome  this  shortage,  the  inference  of  thermospheric  parameters  from  relatively 
plentiful  ionospheric  data  is  aided  by  exploiting  our  knowledge  of  the  intimate  coupling 
between  the  thermosphere  and  ionosphere  using  the  ensemble  Kalman  filtering  (EnKF) 
data  assimilation  system,  constructed  with  the  NCAR  Data  Assimilation  Research 
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Testbed  [Anderson  et  al.,  2009]  and  the  Thermosphere-Ionosphere-Electrodynamics 
General  Circulation  Model  (TIEGCM)  [Richmond  et  ah,  1992].  The  novelty  of  this 
approach  [Matsuo  et  ah,  2013]  allows  taking  advantage  of  the  expansive  increase  in 
recent  years  of  the  ability  of  GPS  systems  to  globally  monitor  ionospheric  parameters. 
Inferred  thennospheric  parameters  are  compared  to  independent  observations  obtained 
from  the  accelerometer  on  board  the  Challenging  Minisatellite  Payload  (CHAMP) 
[Reigber  et  al.,  2002]. 


Finding  1 

“Assimilation  of  COSMIC/FORMOSAT-3  electron  density  data  into  a  coupled 
thermosphere-ionosphere  general  circulation  model  can  extend  the  predictability  of 
thermospheric  mass  density  more  than  3  days” 

The  observing  system  simulation  experiments  (OSSEs)  are  conducted  for  a 
geomagnetically  quiet  12-hour  period  from  0000  UT  through  1200  UT  of  March  20, 
2009,  with  a  one-hour  assimilation  window  and  90  ensemble  members.  Note  that  90- 
member  ensemble  simulations  are  not  centered  at  the  truth  from  which  synthetic 
observations  for  OSSEs  are  sampled,  so  the  model  ensemble  is  not  entirely  free  of  model 
biases.  This  12-hour  data  assimilation  period  is  followed  by  the  72-hour  ensemble 
forecast  period.  Ensemble  forecast  simulations  are  initialized  by  the  assimilation  analyses 
obtained  by  the  OSSEs.  To  examine  the  impact  of  data  assimilation,  additional  ensemble 
simulations  are  run  without  data  assimilation  and  referred  to  as  the  control.  The 
experiment  setup  is  described  in  detail  in  Hsu  et  al.  [2016],  Figure  1  compares  the  neutral 
mass  density  errors  from  three  90-member  ensemble  forecast  experiments,  which  are 
initialized  with  results  from  different  OSSEs,  to  the  control  case.  These  OSSEs  are 
executed  with  different  filters  constructed  with  different  sets  of  the  TIEGCM  state 
variables  being  included  in  the  ensemble  Kalman  filtering  state  vector.  Suppose  f 
represents  a  vector  of  a  given  TIEGCM  model  state  variable  defined  on  the  model  grid. 
The  EnKF  state  vector  for  each  filter  is  composed  of  (a)  the  electron  density  fe-,  the 
atomic  oxygen  ionf0+,  and  the  temperature  fT;  (b)  fe-,  f0+,  fT,  the  atomic  oxygen 
mixing  ratio  f0,  and  the  molecular  oxygen  mixing  ratio  f02;  and  (c)  fe-,  f0+,  fT,  f0,  fo2, 
as  well  as  the  horizontal  neutral  winds  fu,  fv.  Note  that  the  sum  of  the  mixing  ratios  of 
atomic  oxygen,  molecular  oxygen,  and  molecular  nitrogen  is  set  to  1  in  the  TIEGCM,  so 
the  molecular  nitrogen  mixing  ratio  is  adjusted  according  to  the  changes  in  the  atomic 
oxygen  mixing  ratio  f0  and  the  molecular  oxygen  mixing  ratio  f02,  even  thought  it  is  not 
included  in  the  state  vector. 

The  comparison  suggests  that  thermospheric  states,  particularly  temperature,  can  be  well 
inferred  by  assimilating  COSMIC/FORMOSAT-3  data  into  the  TIEGCM,  which  leads  to 
an  improvement  of  the  neutral  mass  density  forecast  by  more  than  3  days.  Effects  of  the 
neutral  composition  initialization,  though  smaller  than  the  impact  of  neutral  temperature 
initialization,  are  long-lasting  when  combined  with  the  neutral  wind  initialization. 
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Figure  1.  The  root-mean-squared  error  (RMSE)  of  the  neutral  mass  density  over  the 
course  of  the  72-hour  ensemble  forecast  period.  RMSE  of  the  control  experiment,  in 
which  no  data  is  assimilated  during  the  preceding  12-hour  assimilation,  is  shown  in  gray. 
RMSE  of  the  90-member  ensemble  forecast  experiment  is  shown  in  light  green,  yellow, 
and  dark  green,  respectively,  for  model  initialization  with  the  three  different  filters 
described  in  the  main  text. 


Finding  2 

“Thermospheric  mass  density  can  be  estimated  from  COSMIC/FORMOSAT-3 
electron  density  data” 

Actual  COSMIC/FORMOSAT-3  electron  density  profile  data  are  assimilated  into  the 
TIEGCM  for  a  geomagnetically  quiet  period  of  23-25  June,  2008  using  the  ensemble  data 
assimilation  experiment  setup  described  in  Lee  et  al.  [2012].  In  this  experiment,  the 
EnKF  state  vector  is  composed  of  fe-,  fT,  fu<  fv<  f o>  fo2  •  A  one-hour  assimilation  window 
and  90  ensemble  members  are  used.  To  validate  the  neutral  mass  density  analysis 
resulting  from  this  assimilation  experiment,  the  result  is  compared  to  independent 
CHAMP  in-situ  mass  density  observations.  During  this  period,  the  CHAMP  satellite 
orbits  at  330-360  km  along  the  noon-midnight  local  time  plane  30  times. 
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Figure  2.  The  root-mean-squared  difference  (RMSD)  of  CHAMP  observations  with  the 
neutral  mass  density  analyses  and  the  control  simulation  is  shown  in  red  and  blue, 
respectively,  with  the  noon  and  midnight  orbit  displayed  by  solid  and  dashed  line. 


'/ 
'  / 


- Control/Midnight 

- Assimilation/Midnight 

- Control/Noon 

- Assimilation/Noon 


0.8 


1.2 

LKg/m3] 


1.4 


1.6 
x  10" 


The  squared  difference  between  the  CHAMP  neutral  mass  density  and  the  density 
analysis  evaluated  at  the  CHAMP  location  is  computed  and  sorted  for  the  noon  and 
midnight  orbits  with  respect  to  10°  latitude  bins,  and  summarized  as  the  RMSD  as  shown 
in  Figure  3.  For  the  noon  orbit,  the  neutral  mass  density  inferred  from  the  assimilation  of 
COSMIC/FORMOSAT-3  electron  density  observation  agrees  better  with  the  CHAMP 
observations  than  the  ensemble  control  simulations.  For  the  midnight  orbit,  the  neutral 
mass  density  specification  is  not  improved  by  data  assimilation  over  the  control 
simulation. 


Finding  3 

“Thermospheric  composition  can  be  estimated  from  GPS-based  electron  density 
data” 

Predictability  of  the  ionosphere  can  be  extended  considerably  by  estimating  neutral 
compositions  from  GPS-based  ionospheric  observations  by  the  coupled  thermosphere- 
ionosphere  data  assimilation  approach  developed  in  this  project.  The  fact  that  the 
ionospheric  forecast  is  aided  by  the  inference  of  thermospheric  compositions  is 
confirmed  within  the  context  of  OSSEs  for  the  COSMIC  I  mission  [Hsu  et  al.,  2014; 
2016]  and  for  the  COSMIC  II  mission  [Lee  et  al.,  2013].  Furthermore,  Chartier  et  al. 
[2016]  and  Chen  et  al.  [2016]  adopted  the  approach  to  assimilation  of  actual  ground- 
based  TEC  data  during  storm  periods,  and  showed  that  the  approach  indeed  extended  the 
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ionospheric  forecast  capability  of  a  coupled  thennosphere-ionosphere  general  circulation 
model.  Though  these  are  indirect  evidences,  thennospheric  compositions  can  be 
estimated  from  both  space-based  and  ground-based  GPS  observations  to  the  extent  that 
results  in  a  considerable  improvement  of  the  ionospheric  forecast  by  up  to  24  hours. 
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